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ABSTRACT 

We demonstrate here using both single-dish and interferometric observations that CO 
molecules are an excellent kinematic tracer, even in high-mass galaxies, allowing us to in- 
vestigate for the first time the CO Tully-Fisher relation of early-type galaxies. We compare the 
Tully-Fisher relations produced using both single-dish and interferometric data and various 
inclination estimation methods, and evaluate the use of the velocity profile shape as a crite- 
rion for selecting galaxies in which the molecular gas extends beyond the peak of the rotation 
curve. We show that the gradient and zero-point of the best-fit relations are robust, indepen- 
dent of the velocity measure and inclination used, and agree with those of relations derived 
using stellar kinematics. We also show that the early-type CO Tully-Fisher (CO-TF) relation 
is offset from the CO-TF of spirals by 0.98 ± 0.22 magnitude at if s -band, in line with other 
results. The intrinsic scatter of the relation is found to be «0.4 magnitudes, similar to the level 
found in the spiral galaxy population. Next generation facilities such as the Large Millimeter 
Telescope (LMT) and the Atacama Large Millimeter/Sub-millimeter Array (ALMA) should 
allow this technique to be used in higher-redshift systems, providing a simple new tool to trace 
the mass-to-light ratio evolution of the most massive galaxies over cosmic time. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: kinematics 
and dynamics - ISM: kinematics and dynamics - galaxies: spiral - galaxies: structure 
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1 Introduction 

The Tully-Fisher relation (TFR; Tully & Fisher 1977 1 of spiral 
galaxies has proved itself to be one of the most important correla- 
tions in extragalactic astrophysics. For example, its use as a distance 
measure is vital in extending the cosmic distance ladder, allowing 
the scale of structures in the nearby universe to be determined and 
studied. The underlying cause of this relation between luminosity 
and rotational velocity is usually interpreted as the product of a rel- 
atively constant total (luminous plus dark) mass-to-light ratio (M/L) 
in the local spiral galaxy population (Gavazzi 1993, Zwaan et al. 
1995 ), and hence a strong coupling between dark and luminous 
mass. Studying the slope and zero-point of the TFR is thus also 
a powerful probe of the M/L evolution of galaxies (e.g. [Phillipps 
[T^89"l|Sprayberry et al.|1995||Bell & de Jong|2001). 

The TFR as introduced by |Tully & Fisher| ( |1977] > uses rota- 
tion velocities derived from HI line-widths in spiral galaxies, which 
should approximate the true projected circular velocity of the galax- 
ies as long as the HI distributions are relaxed and reach into the flat 
parts of the galaxy rotation curves. Over the years, however, the 
TFR has gradually been recognised as denoting the empirical re- 
lationship between the luminosity and rotation velocity of galaxies 
(generally disc galaxies) as measured with a variety of kinematic 
tracers, using differing techniques at various wavelengths. 

Studying the TFR in lenticular and elliptical galaxies (col- 
lectively referred to as early-type galaxies; ETGs) is problematic, 
as they do not all possess extended relaxed atomic gas distribu- 
tions. Indeed, it has been shown that even where HI is present in 
early-type galaxies it can be disturbed, and hence the measured 
line-widths are in some cases unrelated to galaxy properties (e.g. 
|Morganti et al.|2006||Williams et al.|2010| l. Although widespread, 
the ionised gas is generally faint, with significant pressure support; 
and is thus not ideal either (Bertola et al. 1995). The TFR of early- 
type galaxies is nevertheless important, as these galaxies are be- 
lieved to have turbulent formation histories and the TFR may give 
clues about their assembly and evolution. For example, many au- 
thors have suggested that SO galaxies have avoided violent inter- 
actions, and are the faded descendants of high-redshift spirals (e.g. 
|Dresslerp9 80 ; Dressi er" et al.|199"7) . In this scenario, SO galaxies 
become dimmer whilst keeping the same dynamical mass, leading 
to an offset TFR. 

Due to the problems listed above, stellar tracers of galactic 
rotation are generally used in ETGs. Due to the importance of pres- 
sure support in these systems, however, stellar dynamical model- 
ing or an asymmetric drift correction is required to extract the true 
circular velocities. This adds additional systematic uncertainties to 
the already challenging stellar kinematic observations (which must 
reach sufficiently large radii). Recent analyses of ETGs have sug- 
gested that SO galaxies do indeed have a measurable offset from 



the spiral TFR, of around 0.5 - 1.0 mag at _K-band (Neistein et al. 


1999 Bedregal et al. 2006 ; Williams et al. 


2010). Magorrian & Bal- 


lantyneH2001J,|Gerhard et al.|p001) and 


Rijcke etal.H2007|l have 



all considered extending this approach to construct TFRs that in- 
clude elliptical galaxies, and find that these are also offset from the 
spiral TFR, by 0.5 to 1.5 mag at optical wavelengths. 

In this work we consider the use of carbon monoxide (CO) as 
a tracer of the circular velocity, of fast-rotating early-type galaxies. 
CO is thought to be free of many of the problems that beset HI in 
early-types. As part of the ATLAS 3D survey of ETGs i Cappellari 



hereafter Paper IV) , with ~22% of early-type galaxies in the lo- 
cal volume (<42 Mpc distant) containing a substantial molecular 
gas reservoir. The detection rate of molecular gas is also indepen- 
dent of galaxy luminosity and mass (the two most important TFR 
parameters), providing a direct and unbiased probe of the poten- 
tial in high-mass galaxies, which are often HI poor (Serra et al., 
in preparation). The molecular gas is likely to be relaxed in most 
galaxies under study, even in clusters where HI may be undergoing 
ram-pressure stripping ( jToribio & Solanes 2009), due to the short 
dynamical timescales in the central parts of the galaxies, where it is 
usually fou nd (e.g. |Wrobel & Kenney|1992 |Young|2002| |Young| 
[eraLl[2008l |Crocker et al.||2008||2009||2010| >. Therefore CO pro- 
vides a powerful, directly observable measure of the circular veloc- 
ity of galaxies of all masses and morphological types, irrespective 
of enviroment. Additionally, the small beamsizes of the mm-wave 
single-dish telescopes used to detect CO eliminate source confu- 
sion, at least in nearby galaxies, while current interferometers rou- 
tinely yield arcsecond angular resolutions. 

The possibility of using CO line-widths to investigate the 
Tully-Fisher relation was first explored by Dickey & Kazes ^1992) , 
and this method has since been used by various authors to investi- 
gate the CO TFR of spiral and irregular ga laxies (e.g. [ Sch oniger 
& Sofue||1994)|Schoeniger & Sofue|[l997l |Tutui & Sofue||1997| 



|et al.|[20TT| hereafter Paper I) we have found that molecular gas 
is reasonably abundant in early-type galaxies (Youn g et al.||201l| 



Lavezzi & D ickey 1998 , Tu tui et al.|2001| > and quasars (Ho 2007). 
In spiral galaxies, the CO velocity widths obtained are directly com- 
parable to those found in HI (Lavezzi 1997). [Young et ak| ( |2008} 
have also demonstrated that CO is a good tracer of the circular ve- 
locity in a small sample of ETGs, from a comparison with detailed 
dynamical modeling of the stellar kinematics, but this has yet to be 
demonstrated for a larger sample. 

One possible complication introduced when using CO as a dy- 
namical tracer is that the molecular gas in early-types is usually 
confined to the inner regions. This means that we are only able to 
probe the kinematics where dark matter does not yet play a signif- 
icant role. This is the opposite of HI TF analyses, which are gen- 
erally probing HI at large radii, in the dark matter-dominated part 
of the rotation curve. As such, a TFR measured in the central re- 
gion could be considered more similar to the Faber-Jackson (central 
velocity dispersion - luminosity) relation JFaber & Jackson 1976} . 
The ETG CO TFR will therefore only be directly comparable to 
TF results from HI if the 'disc-halo conspiracy' , where dark matter 
flattens the rotation curve at a similar velocity to that found (from 
luminous material) in the inner regions, also holds for early-type 
galaxies ( |Kent|1987[[Sancisi|2004||Gavazzi et al.|2007| >. 

Some authors have suggested that there is a change in the slope 
of the TFR for high-mass disc galaxies, brighter than an absolute 
if-band magnitude of »-23.75 (Peletier & Willner|1993[ |Verhei-| 
pn]|2001| |Noordermeer & Verheijen||2007) . A break in the TFR 
at a similar position is also found for elliptical galaxies ( |Gerhard| 
|et al.|200l[|Rijcke et al.|2007) It has been suggested that this break 
occurs because many massive galaxies have declining circular ve- 
locity profiles, whereas low mass galaxies have relatively flat cir- 
cular velocity curves. In other words, the aforementioned disc-halo 
conspiracy is not perfect in ETGs, and many have a local rotation 
velocity peak at small radii. The radius where one measures the 
velocity hence becomes important (see Noordermeer & Verheijen] 
|2007} . Measures of maximal rotation produce a much larger break 
than measures of the asymptotic rotation velocity, supporting this 
interpretation. Clearly, the existence of a break in the TFR could 
lead to systematic biases when deriving distances or probing galaxy 
evolution. Fortunately, however, we show in this paper that our CO 
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velocities are consistent with measures of the rotation beyond the 
peak of the galaxy circular velocity curves. 

This work, to the best of the author's knowledge, represents 
the first attempt to create a CO TFR for early-type galaxies. One 
can identify two major pitfalls that need to be overcome for this 
approach to be successful. Firstly, as discussed above, it has been 
shown that the CO in ETGs is often very centrally concentrated, and 
hence in some galaxies it may not reach beyond the peak of the rota- 
tion curve. Identifying such cases is critical to obtain a useful TFR. 
Secondly, it has been shown that the molecular gas in early-type 



galaxies is often misaligned with respect to the stars (e.g. [Young 



2002p chinnerer & Scoville 2002| |Young et al.|2008[|Crocker et al. 
2009 a; |Crocker et al.|2008p ). In these cases the inclination of the 
stars is not useful, and we require an estimate of the inclination of 
the molecular gas itself in order to de-project its observed rotation 
velocity. 

Our goal in this paper is therefore to explore and demonstrate 
that, despite the potential pitfalls highlighted above, ETGs do ap- 
pear to follow a robust luminosity-rotational velocity relation, con- 
sistent with that measured using other tracers at the same radii. Al- 
though technical aspects will certainly be improved and the inter- 
pretation of this relation remains uncertain, an ETG CO TFR is a 
tantalising and promising tool for galaxy evolution studies, worthy 
of further consideration. In Section 2 of this paper we discuss the 
data used in this work, outlining the observations in Sections 2.3 
and 2.4, the methods for estimating the inclination in Section 2.5, 
and our method for extracting the circular velocity from observed 
CO line-widths in Section 2.6. In Section 3 we compare these ve- 
locities with other measures such as circular velocities derived from 
dynamical models, and show that it is possible to use simple cri- 
teria to select galaxies where CO is a good tracer of the circular 
velocity beyond the peak of the galaxy rotation curve. In Section 4 
we present our ETG CO TFRs, and explore how different data and 
inclination estimates affect the resulting relations. In Section 5 we 
more fully discuss our results and compare them with previous TFR 
results derived using other tracers. We summarise our conclusions 
in Section 6. 

2 Data 
2.1 Sample 

The sample used here is composed of the galaxies detected in CO(l- 

0) emission in Paper IV. CO(l-O) and (2-1) were observed in all 
galaxies from the complete, volume-limited ATLAS 3D sample of 
260 lenticular and elliptical galaxies (see Paper I for full details). 
Upper limits for non-detections are typically in the range 10 7 -10 8 
Mq, depending on the distance, assuming a standard CO-to-H2 
conversion factor. The detected galaxies used here range in molecu- 
lar gas masses from 10 7 to 10 9 Mq, with molecular gas mass frac- 
tions between 10 -4 and 10 _1 Mq/Lkb, where Lk s is the AVband 
luminosity from the Two Micron All Sky Survey (2MASS) l |Jarrett| 
et al. 2000; Skrutskie et al. 20061. Only fast-rotating galaxies were 
detected. 

Of the 52 detections by Paper IV, 4 are only detected in CO(2- 

1) , and a further 8 had an insufficient signal to noise ratio to per- 
form the analysis required here, leaving 40. In addition to these we 
include the ETGs IC2099, NGC4292 and NGC4309, which were 
detected in CO as part of the initial ATLAS 3D survey, but were later 
removed from the sample because they are too faint at if s -band to 
meet the final sample selection criteria. This should not affect the 
TFR. This leaves a total of 43 galaxies, the properties of which are 
listed in Table Q] 



2.2 Photometric data 

The A'-band magnitudes used in this paper are A^totai from 
2MASS ( |Jarrett et al.(2000l|Skrutskie et al.|2006) . These A" 8jt otai 
magnitudes are measured over large apertures, to include the to- 
tal flux from the galaxy using the techniques developed in Kron 
( |1980| l and curves-of-growth (see Jarrett et al. 2000 for further de- 
tails). They have been widely used in the astronomical community 
and are found to be robust to ~0.1 mag (Noordermeer & V erheijeri| 
[2007} . 

2.3 CO single-dish data 

As reported in Paper IV, the Institut de Radioastronomie Mil- 
limetrique (IRAM) 30m telescope at Pico Veleta, Spain, was used 
for simultaneous observations of CO(1-0) and CO(2-l) in our 
galaxies. The primary beam FWHM is 23" and 12" for CO(1-0) 
and CO(2-l) respectively. The filterbank back-end gave an effec- 
tive total bandwidth of 512 MHz (~1330 km s _1 ) and a raw spec- 
tral resolution of 1 MHz (2.6 km s _1 ) for CO(1-0). The system 
temperatures ranged from 190 to 420 K for CO(1-0). The time on 
source was interactively adjusted so that the final, co-added CO(l- 
0) spectrum for each galaxy had a rms noise level of ~3.0 mK TJ 
(~19 mjy) per binned 31 km s" 1 channel. For further details, see 
Paper IV. 

The CO(1-0) and CO(2-l) spectra for each detected galaxy, 
binned to 31 km s -1 channels, were analyzed to find the velocities 
at which the flux drops to 20% of the peak height. Each spectrum 
was read into a routine which, working outwards from the galaxy 
systemic velocity, locates the first channel where the velocity drops 
to 20% of the peak value. The measurements were inspected by eye 
to ensure noise peaks and troughs were not affecting the result. This 
is analogous to the method used by Tully & Fisher ( 1977 ) for mea- 
suring HI line-widths. The resulting velocity width is henceforth 
denoted W2o- If the molecular gas is a good tracer of the circu- 
lar velocity in our galaxies, then the line-width should approximate 
twice the projected rotational velocity, if the molecular gas distribu- 
tion reaches beyond the peak of the galaxy rotation curve ( |Dickey| 
|&Kazes|1992} . 

Figure [T] shows a comparison between W20 line- widths mea- 
sured from CO(1-0) and CO(2-l) spectra. One might expect it to 
be easier to measure line-widths using CO(2-l), as Paper IV have 
shown that this transition is usually 1^1 times brighter than CO(l- 
0) (due to beam dilution and/or intrinstic effects), but Figure [T] 
shows that the agreement between the two measurements is gen- 
erally good. However, almost all the scatter is below the 1 : 1 line. In 
these cases, the molecular gas likely extends further than the «12" 
FWHM of the CO(2- 1 ) beam, and the CO(2- 1 ) line- width measure- 
ments are systematically biased low. In Section pi. 1.3 [ we in fact 
show that in some cases even the «23" CO(1-0) beam is not ex- 
tended enough to retrieve the full velocity width. The number of 
galaxies for which this is a concern is small, however, and in the 
rest of this paper we will therefore use line-widths measured from 
the CO(1-0) spectra. 

The instrumental dispersion is small compared to our 31 km 
s _1 channels, so we did not correct the line-widths for it. We also 
chose not to use corrections for turbulence or line-broadening, as 
these are generally derived by comparison with HI line-widths in 
spirals ([Lavezzi & Dickey[[T^98l [Tully & Fouque||1985| [Tutui &| 
|Sofue|1999 1 and hence may not be applicable to early-type galax- 
ies. The velocity dispersion in the gas is expected to be small (e.g. 
Okuda et al. 2005 1. In Section [3~7] we attempt to quantify any biases 
introduced by ignoring these corrections. The derived W20 velocity 
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Figure 1. A comparison between CO(l-O) and CO(2-l) single-dish line- 
widths, measured as described in Section |2~3"1 The solid line shows the one- 
to-one relation. Linewidths are discretized due to the 30 km s" 1 channel 
width. 

widths are estimated to be robust to within half a channel width, 
~15 km s" 1 , and are listed in Table [T] 

2.4 Interferometric data 

As part of the ATLAS 3D survey, all CO(l-O) detections with an 
integrated flux greater than about 19 Jy km s _1 that do not have 
interferometric maps available in the literature are to be observed 
with the Combined Array for Research in mm-wave Astronomy 
(C ARMA; [Bock et al.|2006) . Full details of this ongoing interfero- 
metric survey can be found in Alatalo et al. (in preparation), but we 
summarize the observations here. 

Observations of the sample galaxies have been ongoing since 
early 2008, mainly in the D-array configuration, providing a spatial 
resolution of 4-5". CO(l-O) has so far been observed using narrow- 
band correlator configurations, providing at least 3 raw channels 
per 10 km s _1 binned channel whilst ensuring adequate velocity 
coverage for all galaxies. Bright quasars were used to calibrate the 
antenna-based gains and for passband calibration. The data were 
calibrated and imaged using the 'Multichannel Image Reconstruc- 
tion, Image Analysis and Display' (MIRIAD) software package 
^Sault et a1"1|1995fr . Total fluxes were compared with the IRAM 
30m single-dish observations to ensure that large proportions of the 
fluxes were not being resolved out. 

A total of 22 galaxies included in this work have been ob- 
served with CARMA so far (IC 0676, IC 719, IC 1024, NGC 1222, 
NGC1266, NGC 2764, NGC 2824, NGC 3626, NGC 3665, 
NGC4119, NGC4292, NGC4324, NGC4429, NGC4435, 
NGC4694, NGC4710, NGC5379, NGC6014, NGC7465, 
PGC058114, UGC06176 and UGC09519). We also include 
here the galaxies for which data is already available from the 
literature, mostly from SAURON survey ( |de Zeeuw et al.|[2002) 
follow-ups. These are NGC 3032, NGC 4150, NGC 4459 and 
NGC 4526 \ Young et al.|2008); NG C 2685 ( |Schinnerer & Scoville| 
[2002} ; NGC 2768 ( |Crocker et aL]|2008l >; NGC 0524, NGC 3489, 
and NGC 4477 ( |Crocker et al.|2010| >. This makes for a total of 31 
sample galaxies that have interferometric data. 

The data cubes were summed spatially to measure total fluxes 
and revised values for W20, serving as consistency checks on the 
values derived from the single-dish data. The large primary beam 
of the arrays also allow us to make better measurements for the 



small subset of galaxies in which a substantial part of the CO dis- 
tribution was missed by the 23" beam of the IRAM 30m telescope, 
such as in IC676, NGC 4324, NGC 4477, NGC 47 10, NGC 7465 
and PGC058114. This can be due to pointing errors or molecu- 
lar gas distributions that extend beyond the beam. The integration 
time for the single-dish CO observations was adjusted to obtain a 
fixed noise level- and as such galaxies which are only detected at 
«5er will have rather uncertain W20 line-widths. The interferomet- 
ric observations in these cases (NGC 2685, NGC 2768, NGC 3489, 
NGC 4477) will provide a better constraint on the line-widths. The 
revised W20 values are listed in Table [T] but for consistency these 
velocities are only used in the TFRs presented in Sections |4. 1.3| and 

MM 

2.5 Inclination correction 

The measured quantity W20 is a projection of the gas velocity into 
the line-of-sight. If this is to be used for a TF analysis it must be 
deprojected. Many methods for inclination measurement are avail- 
able. Starting with the simplest, we compare various method below 
in order of increasing complexity, to allow future CO TF surveys 
to select the optimum method for their needs. The methods pre- 
sented in Sections |2.5. 1 1 and |2.5.2| use only single-dish molecular 
data, whilst the method presented in Section [2.5.3| requires interfer- 
ometric maps. 

2.5.1 Galaxy axial ratio 

A rough measurement of inclination can be obtained by estimating 
the axial ratio of the stellar distribution of the host galaxy from 
imaging data: 




where q is the ratio of the semi-minor (6) to the semi-major (a) axis 
of the galaxy, and go is the intrinsic axial ratio when the galaxy 
is seen edge-on (qo = c/a). qo is often assumed to be 0.2 in 
disc galaxies JTully & Fisher|1977[ >, but early-types can have large 
bulge-to-disc ratios leading to a large uncertainty in any assumed 
value of go ■ 

Various lines of enquiry suggest a mean qo value of ~0.34 for 
the fast-rotators in the ATLAS 3D sample (Weijmans et al., in prepa- 
ration), which we adopt here. Intrinsic scatter around this value will 
introduce an artificial increase in the TFR scatter, but the effect is 
very small. Indeed, we stress that this go-related inclination correc- 
tion is only significant in highly inclined galaxies, where fortunately 
inclination errors have the smallest effect on the de-projected veloc- 
ities and thus the derived TFR. For example, assuming go = for 
the current sample would only change the zero-point of the derived 
TFR by 0.02 dex along the velocity axis, and it does not change the 
gradient. 

It is important to note that using this method to correct the CO 
velocity widths for a TF analysis implicitly assumes that the CO is 
distributed in the same plane as the galaxy major-axis. Recent work 
by Davis et al., (in preparation) however suggests that a sizable pro- 
portion of early-type galaxies have molecular gas misaligned with 
respect to the stellar kinematic axis, so in these cases one expects 
this method to introduce additional artificial scatter in the TFR. 

The value of ib/ a obtained by assuming go=0.34 for each 
galaxy is listed in Table[T| These values were calculated using Equa- 
tion[T|with the mean galaxy axial ratio reported in the NASA/IPAC 
Extragalactic Database (NED). These values are the mean of the 
measurements from the Sloan Digital Sky Survey (SDSS) in the r 
band, 2MASS at K s band < |Skrutskie et al.|2006) , and blue and red 
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filters from both the Third Reference Catalogue of Bright Galaxies 
(RC3; [de Vaucouleurs et al.|1991|, and the Uppsala General Cata- 
logue of galaxies (UGC/POSS; |Nilson|1973"l > where available. The 
exceptions to this are NGC 3665, where the axial ratio is based only 
on the superior SDSS imaging, and PGC058114, where the mean 
value of the axial ratio from the HyperLeda database ( Patu rel et al.| 
20031 has been used, due to a larger number of available measure- 
ments. When these inclinations are compared to those estimated 
from the axial ratios of Krajnovic et al. (20101 (Paper II; derived by 
calculating the moments of inertia of the surface brightness distri- 
bution from the SDSS and INT r-band images) for the galaxies in 
the ATLAS 3D sample they are found to agree well. 

2.5.2 Dust axial ratio 

It should be possible to obtain more accurate inclination estimates 
from fitting ellipses to dust highlighted in unsharp-masked optical 
images of the galaxies. This method has several advantages over 
galaxy axial ratios. Dust distributions typically have very small ver- 
tical scale-heights, hence Equation[T]wifh qo=0 should yield a good 
estimate of their inclinations. It has also been shown that CO and 
dust in early-type galaxies are usually spatially coincident ( |Young| 
|et al.|2008| [Crocker et al.|2008||2009|[20T0l >, hence the inclination 
of the dust (idust) should trace the true inclination of the molecular 
gas. This is especially useful where the molecular gas and stars are 
misaligned. 

Unsharp-masked dust maps were created for all the galaxies 
in this work, preferentially from archival Hubble Space Telescope 
(HST) images, or where these were n ot available from SDSS g- 
band |Adelman-McCarthy et al.|2 008 1 or ATLAS 3 D Isaac Newton 
Telescope Wide Field Camera (INT-WFC) r-band images (Scott et 
al., in preparation). Ellipses were then fitted by eye to the resulting 
maps and the inclination of the dust calculated from Equation[T]with 
go=0. Four example galaxies are shown in Figure|2] 

One of the galaxies (NGC 4753) is well studied in the liter- 
ature, showing many thin filamentary dust features that are well 
fitted by a model featuring an inclined disc, twisted by differential 
precession, with an average id us t=75 ± 3° ( Stei man-Cameron et al.| 
19921. We adopt this value here. A full list of idust values can be 
found in Table Q] 

2.5.3 Molecular gas modeling 

Tilted-ring analyses have been used to determine the geometric 
and kinematic parameters of neutral and molecular gas discs in 
interferometric data for many years (e.g. |Rogstad et al.||1974| 
IChristodoulou et al.|1988l|Koribalski et al.|1993) . In this work, we 
have made use of the Groningen Image Processing System (GIPSY) 
add-on package TiRiFiC, described in |J6zsa et al.| (|2007 1, to fit 
tilted-ring models to observed data cubes. Due to the small num- 
ber of independent synthesized beams across our typical sources, 
an unwarped disc model was fitted where the inclination, position 
angle, systemic velocity and kinematic centre of the gas were varied 
globally, and only the velocity and surface brightness were allowed 
to vary radially (i.e. for each ring). Full details of the fits are re- 
ported in Davis et al., (in preparation) but the inclination measures 
(i mo i) used in this work are listed in Table[T] The error quoted for 
each galaxy is the maximum inclination difference found if each 
side of the velocity distribution is supplied as an input to the fitting 
program separately. 

We have attempted the analysis described above for all the 
mapped galaxies in our sample where with have adequate spatial 
resolution, with one exception where we felt the literature value 
was sufficient. NGC 2685 has dust lanes which look polar on the 




(c) NGC 3626 



(d) NGC 5379 



Figure 2. Four examples of unsharp-masked images picking up dusty spiral 
structures and rings. Overlaid in blue are the best fit ellipse from which we 
estimate the dust inclination. 



north-east side, and are coincident with HI, Ha and CO emission 
( jSchinnerer & Sc oville 2002J. This system has been revealed to 
feature a coherent, extremely warped disk (Jozsa et al. 2009), for 
which we adopt the published average inclination value of i mol =69 
±9°. 

2.6 W20 linewidths 

The velocity widths calculated as described in Sections [2.3| and [2~4| 
are an approximation to (twice) the maximum rotation velocity of 
the CO. To be useful, these CO measurements must reach beyond 
the peak of the galaxy rotation curves. Young et al. ( 2008 1 suggest 
that molecular gas in early-types is often centrally concentrated and 
hence this requirement is not trivially fulfilled. 

One sign of a tracer having reached a flat part of a galaxy ro- 
tation curve is a classic 'double-horned' velocity profile. Lavezzi & 
|Dickey| ( |l997} discuss in detail the problem of recovering velocity 
widths from CO observations, given the variety of profile shapes 
typically seen in millimeter observations. They show that CO pro- 
file shapes are affected by the extent of the emitting gas, the beam 
size of the telescope, pointing errors and the optical depth of the 
molecular material. CO line-widths are shown to accurately retrieve 
Vsat in the majority of cases where the flaring parameter R2,s (the 
ratio of the velocity widths measured at 20 and 50% of the peak) is 
less than 1.2 ( |Lavezzi & Dickey|1997l[T99"8") . This criterion effec- 
tively selects galaxies that show a boxcar, or double-horned, pro- 
file shape with sharp edges and rejects those with a more Gaussian 
profile. We follow this methodology in this work, exclusively using 
galaxies with double-horned/boxcar/sharp edged profiles, which we 
call 'boxy' from here on in, assuming that in these systems the mea- 
sured gas velocity width will approximate the velocity beyond the 
peak of the circular velocity curve. 

The limited signal-to-noise ratio of the spectra we possess, 
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however, means that the value of the flaring parameter discussed 
above is often rather uncertain. We performed several other tests, at- 
tempting to find a good automatic way of classifying boxy velocity 
profiles. These included fitting single gaussians and selecting those 
galaxies that shown significant structure in the residuals, and fitting 
Gauss-Hermite polynomials and looking for a significant kurtosis 
(I14 < 0). Both of these methods worked well for wide lines and 
high signal to noise cases, but in the case of narrow lines, or low 
signal to noise the fits were not well constrained. In these cases the 
different methods often produced contradictory results, and/or did 
not agree with a careful classification by eye. In future surveys with 
better data quality, or large numbers of objects (so one can afford 
to remove low signal to noise detections) these methods will likely 
be highly useful, but in this current work, with a limited number of 
objects we prefer a careful classification by eye. 

Even by eye, some cases are hard to classify, especially when 
the profiles are narrow as the galaxies are fairly face on, or if the 
signal-to-noise ratio is low. We cautiously use only those galaxies 
with clear boxy profiles here, flagging the 12 galaxies with uncer- 
tain profiles shapes. 

Using this profile shape criterion is likely to introduce some 
uncertainties, for example galaxies with a molecular ring will have 
double-horned profiles without the CO necessarily having reached 
beyond the peak of the rotation curve. Similarly, galaxies that have 
strong resonances may have their gas distribution truncated at the 
turnover radius, causing them to not display a classic boxy shape 
even though the edge of the distribution is rotating at the turnover 
velocity. In Section ITT] we use the dynamical models available to 
us from the ATLAS 3JJ survey to quantify these biases. 

2.6.1 Combining interferometric and single-dish data 

Finally, we can combine our single-dish and interferometric obser- 
vations, using all available data to increase the number of galax- 
ies with high-quality measurements available. When available, we 
thus always use interferometric velocity widths and, when well 
constrained, interferometric inclinations (quantitatively, we select 
galaxies with inclination errors less than 6°). For the remaining 
galaxies we preferentially use dust inclinations, or failing this, bla 
inclinations. We can also use the interferometric data when avail- 
able to help ascertain which galaxies have reached beyond the peak 
of the rotation curve, in conjunction with the usual single-dish boxy 
criterion. This hybrid measure represents our best estimate of the in- 
clination and velocity width. The adopted best inclinations are listed 
in Tabled 

3 Comparison with other velocity measures 
3.1 Model circular velocity curves 

As alluded to in Section [2~6| it has yet to be established for a large 
sample of early-type galaxies that single-dish CO is a good tracer 
of the circular velocity. One would like to test this , and attempt 
to quantify the biases introduced by using the profile shape as a 
proxy for reaching beyond the peak of the circular velocity curve. 
As Williams et al. (2010) discuss, one must be careful when com- 
paring velocity measures derived from different methods, as they 
are likely to have significant systematic differences. However, it is 
important to ascertain that the CO line-widths vary systematically, 
in a similar way to other rotation measures. 

The gravitational potential of ETGs can also be significantly 
different from that of spirals. ETGs often have circular velocity pro- 
files with high peaks at small radii that decline before flattening out. 
The W20 line-width is, in general, more sensitive to the peak ve- 
locity of the molecular gas, rather than the flat part of the rotation 



curve, unless the gas disk is sufficiently extended. This would lead 
to an overestimation of the velocity width. This problem has been 
discussed in detail by Noordermeer & Verheijen ( 2007 1, who find 
that using the HI peak velocity in massive early-type spirals results 
in an offset TFR. They report that the asymptotic velocity is a bet- 
ter measure of the total potential, and results in spirals of all masses 
lying on a single TFR. 

An estimate of the amount by which we may be overestimat- 
ing the velocity widths can be obtained using the circular velocity 
curves produced by Williams et al. (2009) from stellar dynamical 
modeling of their SO galaxies. The median value of the circular ve- 
locity curve after the peak out to the maximum extent of the data, 
where the model is well constrained, is adopted as a measure of the 
velocity of its flat part (in the case of [Williams eTal] p009) l then- 
data extends to ~3 Re). When tested on the SO circular velocity 
curves from Williams et al. (2009), this median performs well in 
picking out the velocity of the flat section. The mean ratio of the 
peak velocity to the flat velocity is 1.15, or 0.06 dex. When trans- 
formed to an offset in luminosity using the parameters of the best-fit 
SO TFR of Willi ams et al.| ( [2"0"T0| >, this effect corresponds to an off- 
set of upto 0.5 mag at _fT s -band. This is approximately half of the 
~ 1 mag offset observed by various authors between the spiral and 
ETGs, and is thus significant if not properly accounted for. 

In order to test the assumption that our ETGs with boxy pro- 
files have gas distributions that reach beyond the peak of the circular 
velocity curve, and to see if our line-widths are overly sensitive to 
the peak velocity, like |Williams7uLlp009l|2TjTo) we utilize the ax- 
isymmetric Jeans anisotropic dynamical modeling (JAM) method 
described in Cappellari (2008). Some examples of the approach, 
using SAURON integral-field kinematics (Emsell em et al.|20 04 1 as 
done here, are presented in |Scott et al.U2009|. For the ATLAS 3D 
survey a Multi-Gaussia n expansion (MGE; |Emsellem et al.|1994) 
was fitted to the SDSS jAbazajian et al.|2009) or INT photometry 
(Scott et al. in preparation). The MGEs were then used to construct 
JAM models for all the 260 ATLAS 3 D galaxies (see |Cappellari et al. 
|2010| a) which were fitted to the SAURON stellar kinematics (Paper 
I). The models have three free parameters, the inclination (jjam), 
the mass-to-light ratio (M/L) assumed to be spatially constant, and 
the anisotropy /3 e = 1 — <J^ja\, which is also assumed to be spa- 
tially constant. From each mass model we have calculated the pre- 
dicted circular velocity curve in the plane of the galaxy. It is worth 
bearing in mind that, where the CO is misaligned from the plane 
of the galaxy, we are thus making implicit assumptions about the 
symmetry of the matter distribution. As these models include no 
dark matter, the circular velocity often declines at large radii, and 
one must be careful to measure the circular velocity at a suitable ra- 
dius (i.e. where we have constraining data) when comparing to the 
observations. 

Analysis of the ATLAS 3D interferometric observations, pre- 
sented in Davis et al., (in preparation), suggests that the average 
radial extent of CO discs in early-type galaxies is around one op- 
tical effective radius (the radius encompassing half of the light), 
hence the circular velocity measured at one effective radius (de- 
noted Vjam.Rc) is a sensible quantity to use for comparison. The 
effective radius used here is a combination of the values from 
2MASS ISkrutskie et~aT1|20T)6) and RC3 ( |de Vaucouleurs et al.| 
[199TJ, as described in Paper I. For galaxies where interferometric 
observations are available, we have also calculated the model cir- 
cular velocity at the maximum extent of the CO, and denote this 
Vjam.Rco • Both values can be found in Table[2] 

We compare the two model circular velocity measures de- 
scribed above with the measured CO line-widths in Figure [3] The 
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Table 1. Parameters of the ATLAS 3 early-type CO-rich galaxies used in this paper. 



Galaxy 


T-type 


W 2 0,SD 


W20,intcr 


*b/a 


^dust 


^tilted— ring 


^bcst 




Dhorn 


Dhorn 


Tel. 


Ref. 






(kins" 1 ) 


(kms" 1 ) 


(deg) 


(deg) 


(deg) 


(deg) 


(mag) 


(SD) 


(best) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


IC0676 


-1.3 


93 


170 


58 ± 3 


- 


69 ± 6 


69 ± 6 


-22.20 


- 


- 


C 


1 


IC 0719 


-2.0 


343 


353 


82 ± 1 


79 ± 1 


74 ± 5 


74 ± 5 


-22.70 


X 


X 


C 


1 


IC 1024 


-2.0 


221 


240 


79 ± 1 


72 ± 1 


- 


72 ± 1 


-21.70 


X 


X 


C 


- 


*IC2099 


-1.5 


93 


- 


83 ± 1 


- 


- 


83 ± 1 


-21.40 


9 


? 


- 


- 


NGC0524 


-1.2 


312 


320 


28 ± 8 


19 ±7 


44 ±28 


19 ± 7 


-24.70 


X 


X 


p 


2 


NGC 1222 


-3.0 


156 


210 


41 ±5 


41 ±3 


- 


41 ± 3 


-22.70 


? 


- 


C 


1 


NGC 1266 


-2.1 


187 


180 


47 ±4 


26 ± 5 


- 


26 ± 5 


-22.90 


- 


- 


C 


1 


NGC 2685 


-1.0 


156 


220 


70 ± 1 


61 ± 1 


70 ± 10 


61 ± 1 


-22.80 


- 


- 


o 


4 


NGC 2764 


-2.0 


312 


310 


65 ±2 


65 ± 1 


76 ± 15 


65 ± 1 


-23.20 


X 


X 


C 


1 


NGC 2768 


-4.4 


322 


360 


36 ± 6 


- 


- 


36 ± 6 


-24.70 


X 


X 


p 


2 


NGC 2824 


-2.0 


312 


310 


56 ± 3 


61 ± 1 


- 


61 ± 1 


-22.90 


X 


X 


C 


1 


NGC 3032 


-1.9 


143 


150 


36 ± 6 


35 ± 3 


46 ± 11 


35 ± 3 


-22.00 


X 


X 


B 


3 


NGC 3182 


0.4 


218 


- 


39 ± 5 


35 ± 3 


- 


35 ± 3 


-23.20 


X 


X 


- 


- 


NGC 3489 


-1.2 


292 


240 


62 ±2 


- 


56 ± 15 


62 ± 2 


-23.00 


7 


X 


P 


2 


NGC 3619 


-0.9 


405 


- 


44 ± 4 


48 ± 2 


- 


48 ± 2 


-23.50 


X 


X 


- 


- 


NGC 3626 


-1.0 


374 


374 


51 ±3 


61 ±2 


67 ± 5 


67 ± 5 


-23.30 


X 


X 


C 


1 


NGC 3665 


-2.1 


624 


630 


90 ± 1 


64 ± 1 


74 ±35 


64 ± 1 


-24.90 


X 


X 


C 


1 


NGC 41 19 


-1.3 


156 


170 


83 ± 1 


67 ± 1 


69 ± 3 


69 ± 3 


-22.60 


? 


- 


C 


1 


NGC 4 150 


-2.1 


234 


238 


50 ± 3 


54 ± 2 


54 ± 3 


54 ± 3 


-21.60 


X 


X 


B 


3 


*NGC4292 


-1.7 


187 


190 


49 ± 3 


50 ± 2 


46 ± 14 


50 ± 2 


-21.54 


? 


X 


C 


1 


*NGC4309 


-1.5 


124 


- 


69 ± 2 


61 ± 1 


- 


61 ± 1 


-20.91 


- 


- 


- 


- 


NGC 4324 


-0.9 


218 


360 


73 ± 1 


64 ± 1 


62 ± 1 


62 ± 1 


-22.60 


7 


X 


C 


1 


NGC 4429 


-1.1 


499 


532 


68 ± 2 


68 ± 1 


60 ± 16 


68 ± 1 


-24.30 


X 


X 


C 


1 


NGC 4435 


-2.1 


405 


380 


48 ±4 


52 ± 2 


- 


52 ± 2 


-23.80 


X 


X 


C 


1 


NGC 4459 


-1.4 


391 


400 


43 ±4 


46 ± 2 


47 ± 2 


47 ± 2 


-23.90 


X 


X 


B 


3 


NGC 4477 


-1.9 


204 


260 


36 ± 6 


26 ± 5 


38 ± 3 


38 ± 3 


-23.70 


7 


X 


P 


2 


NGC 4526 


-1.9 


655 


663 


82 ± 1 


78 ± 1 


64 ± 8 


82 ± 1 


-24.60 


X 


X 


B 


3 


NGC 4684 


-1.2 


249 


- 


77 ± 1 


- 


- 


77 ± 1 


-22.20 


7 


7 


- 


- 


NGC 4694 


-2.0 


93 


65 


69 ± 2 


- 


- 


69 ± 2 


-22.10 


- 


- 


C 


1 


NGC 47 10 


-0.9 


312 


430 


78 ± 1 


88 ± 1 


86 ± 6 


86 ± 6 


-23.50 


X 


X 


C 


1 


NGC 4753 


-1.4 


530 


- 


61 ±2 


75 ± 1 


- 


75 ± 1 


-25.10 


X 


X 


- 


5 


NGC 5 173 


-4.9 


187 


- 


24 ± 9 


- 


- 


24 ± 9 


-22.90 


X 


X 


- 


- 


NGC 5273 


-1.9 


218 


- 


38 ± 5 


38 ± 3 


- 


38 ± 3 


-22.40 


X 


X 


- 


- 


NGC 5379 


-2.0 


124 


120 


80 ± 1 


64 ± 1 


- 


64 ± 1 


-22.10 


- 


- 


C 


1 


NGC 6014 


-1.9 


156 


160 


29 ±7 


22 ± 6 


22 ± 10 


22 ± 6 


-23.00 


X 


X 


C 


1 


NGC 7465 


-1.9 


158 


180 


56 ± 3 


70 ± 1 


58 ± 9 


70 ± 1 


-22.80 


X 




C 


1 


PGC 29321 


0.0 


124 




31 ±7 


38 ± 3 




38 ± 3 


-21.60 










PGC 56772 


-2.0 


249 




64 ± 2 


57 ± 2 




57 ± 2 


-22.00 


7 


7 






PGC 58 114 


-2.0 


202 


240 


71 ± 1 




76 ±30 


71 ± 1 


-21.60 




X 


C 


I 


PGC 61468 


0.0 


218 




51 ± 3 






51 ± 3 


-21.60 


7 


7 






UGC 05408 


-3.3 


156 




31 ±7 






31 ± 7 


-21.90 


X 


X 






UGC 06176 


-2.0 


249 


230 


70 ± 2 


68 ± 1 




68 ± 1 


-22.60 


7 




C 


1 


UGC 095 19 


-1.9 


187 


210 


47 ±4 


41 ±3 


63 ± 7 


41 ± 3 


-22.10 


7 


X 


C 


1 



Notes: Column 1 lists the galaxy name. Galaxies with stars next to their name a re early-types that were observed by Paper IV but are not included in the 
ATLAS 3D sample. Column 2 contains the morphological T-type from HyperLeda jPaturel et al.|2003") . This morphology indicator was not used for the sample 
selection (see Paper I). Columns 3 and 4 list the velocity widths at 20% of the peak flux, derived from single-dish and interferometric data, respectively. 
The error in these quantities was estimated as half the velocity width of an individual channel, 15 km s~ 1 . Column 5 contains the inclination derived from 
galaxy axial ratios. The quoted errors only take into account the error in the measurement of the optical axial ratio, as reported in NED. Column 6 contains 
the inclination derived from ellipse fitting to unsharp-masked dust images. The error is estimated to be 5% on the measurement of both the minor and major 
axes. Column 7 contains the inclinations estimated from tilted-ring fits to the interferometric datacubes, the errors on which are calculated as the maximum 
difference from the best value obtained when using only one half of the galaxy velocity field. Column 8 lists our adopted best inclinations, as described in 
Section [2.6.1| Column 9 contains the iC s -band magnitudes from 2MASS (Skrutskie et al. 2006), c onverted to absol ute magnitudes using th e distance to each 
galaxy adopted for the ATLAS 3D survey in Paper I. These distances are drawn preferentially from |Mei et al.|(2007) and |Tonry et al.|{200l| . We assign to the 
absolute magnitudes an error of ±0.1 mag, taking into account distance uncertainties, which completely overwhelm errors caused by the lack of an internal 
extinction correction. Columns 10 and 11 list the galaxies with boxy CO profiles, as determined from visual inspection of the single-dish spectra (Column 9), 
and a combination of single-dish and interferometric data where available (Column 1 1). A question mark indicates that with the respective data we are unable 
to determine if the profile is boxy, and hence do not include it. The telescope used to obtain the interferometric data is indicated in Column 1 2, where C = 
CARMA, P = PdBI, B = BIMA and O = OVRO. Column 13 lists the relevant references for the interferometric observations, and in two cases (NGC 2685 and 
NGC 4753) references for the inclination measurements. (1) Alatalo et al., (in preparation) (2) Crocker et al. (2010), (3) Young et al. 1 20081, (4) Joz sa et ah] 
(20091 and (5) |Steiman-Cameron et aTHl992) . 
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line-widths of the galaxies with boxy CO profiles correlate well 
with the circular velocities measured at both one effective radius 
(R c ) and the radius of maximum CO extent (Rco)- The offset from 
V jam, Re is 0.01 ± 0.02 dex, the data points displaying an RMS 
scatter of 0. 1 dex. The galaxies that do not have boxy CO profiles 
have line-widths systematically smaller than those predicted from 
the models in most cases, suggesting that their CO has been cor- 
rectly identified as not reaching beyond the peak of the rotation 
curve. The comparison with Vjam.Rco, albeit with smaller number 
statistics, shows that the difference between the model and the CO 
line-widths is again consistent with the models, with a mean off- 
set of 0.002 ± 0.02 dex and a RMS scatter of 0.1 dex. There is no 
systematic behaviour in the residuals for the boxy galaxies in either 
plot. The galaxies that do not have boxy CO profiles are closer to 
the JAM model predictions, but do not lie on the one to one relation 
as expected, with a mean offset of -0.07±0.05 dex. This offset is 
likely because the CO is these objects is very compact, and beam 
effects mean the CO size is overestimated. The JAM model veloc- 
ities have thus been estimated at the wrong radius. The one galaxy 
which does not have a boxy profile, but has a line-width larger than 
the predicted circular velocity, is NGC 1266, recently discovered to 
have a kiloparsec-scale molecular outflow ( Alatal oet al.|20lT"| >. 

To check if our line-widths are overly-sensitive to the peak of 
the rotation curve, we also extract from the JAM models the max- 
imum of the circular velocity curve, denoted V jam, max, and the 
median value of the circular velocity curve, measured between the 
peak and the maximum extent of the SAURON data (~1R C ), denoted 
VjAM,> P eak- As discussed above, defining VjAM,> P eak in this way 
picks out the value of the flat part of the rotation curve well in tests 
based on the data of Williams et al. ( 2009), which extend to ~3 Re. 
Our data is only constrained out to around ~1R C , but the rotation 
curves do not drop wildly in this region. It should however be noted 
that this measure may not be the same as the velocity measured at 
very large radii. The values of V JA M,max and Vj A M,>peak can be 
found in Table|2] 

The top panel of Figure[4]shows a comparison between the ob- 
served CO line-widths, de-projected using dust inclinations (as de- 
scribed in Section [2. 5. 2) , and VjAM.max- The CO line-widths are 
systematically smaller than the maximum of the circular velocity 
curve. The median peak of the JAM circular velocity curves is at 
8.5", well within the single-dish beam, and hence we do not expect 
systematically missed emission to cause this trend. The mean offset 
from the maximum value is -0.03 ± 0.02 dex, broadly similar to 
the offset between the maximum and the flat part of the circular ve- 
locity curve expected from the results of Williams et al. ( 2009}. The 
RMS scatter around this value is 0.09 dex. This suggests that we are 
not overly sensitive to the peak of the circular velocity curve. The 
galaxies which do not have boxy profiles are again offset to much 
lower velocities. The only galaxy whose line-width is inconsistent 
with the maximum value of its predicted circular velocity curve is 
again NGC 1266. 

The bottom panel of Figure[4]compares the observed CO line- 
widths, de-projected using dust inclinations, with VjAM,> P cak- The 
CO measurements for galaxies with double horned profiles are, 
within the errors, consistent with the VjAM,> P oak measure. The 
mean offset from this measure is 0.02 ± 0.02 dex, with an RMS 
deviation of 0.1 dex. The galaxies that were not identified as boxy 
are on average offset from Vj A M,> P cak by -0.1 ±0.04 dex. There 
are several galaxies that do not have a boxy profile, but are con- 
sistent with rotating at Vj A M,> pC ak- These galaxies are NGC 4684, 
PGC 56772 and PGC 61468, all of which have profiles which were 
hard to classify. If interferometric observations were available we 
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Figure 3. Difference between our best CO line-widths de-projected using 
the best inclinations (as described in Section |2.6.1) and twice the JAM 
model circular velocities, measured at one effective radius (top) and the ra- 
dius of maximum CO extent, or 12" whichever is smaller (bottom), plotted 
against the best inclination de-projected CO velocity width. Galaxies with 
boxy profiles are plotted as filled diamonds, while the others are shown as 
open squares. The solid line indicates where the JAM and CO velocities are 
identical. The dashed line and dot-dashed line are fitted offsets from zero, 
for the boxy galaxies and the other galaxies, respectively. The errors on the 
fits are indicated as shaded regions around the best fit lines. 



would hence expect to find that the gas does indeed reach beyond 
the peak of the galaxy rotation curve in these systems. 

Figures [3] and [^demonstrate that, if one selects galaxies with 
boxy profiles, W20 velocity widths provide an good estimate of the 
circular velocity after the inner peak of the rotation curve, whether 
estimated as V jam, Re or VjAM,> P eak- It is worth mentioning that, 
for our sample, a simple cut excluding galaxies with CO linewidths 
< 250km s _1 would have a similar effect as selecting by profile 
shape, removing most of the outliers. Applying an arbitrary cut in 
velocity, the most important TF parameter, could bias the derived 
TFR, especially for a sample selected in a different way, and hence 
we prefer selection by profile shape in this work. 

3.2 HI line-widths 

|Morganti et al.| < |2006| > and |Williams et al.| ( |2010) have shown that 
in some cases HI line-widths are not a good measure of the circular 
velocity for early-type galaxies. We examine this here by retrieving 
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Table 2. Parameters extracted from the JAM model circular velocity curves 
for the ATLAS 30 early-type CO-rich galaxies used in this paper. 



Galaxy 


VjAM,Re 


VjAM,Rco 


VjAM,max 


V J AM, > peak 




(km s 1 ) 


/i — 1 \ 
(km s L ) 


(km s 1 ) 


/i — 1 \ 
(km s 1 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


IC 0676 


117 


116 


118 


116 


IC0719 


232 


232 


233 


218 


IC 1024 


150 


163 


174 


157 


IC 2099 






- 


- 


NGC 0524 


337 


362 


368 


337 


NGC 1222 


177 


177 


236 


177 


NGC 1266 


141 


116 


144 


134 


NGC 2685 


160 




186 


160 


NGC 2764 


199 


190 


207 


196 


NGC 2768 


320 


325 


342 


311 


NGC 2824 


220 


220 


278 


206 


NGC 3032 


129 


158 


238 


135 


NGC 3182 


183 


■ 


198 


178 


NGC 3489 


169 


197 


198 


159 


NGC 3619 


216 


• 


277 


223 


NGC 3626 


173 


187 


256 


174 


NGC 3665 


358 


354 


362 


321 


NGC 4150 


134 


141 


163 


131 


NGC 4270 


185 


- 


193 


182 


NGC 4292 




- 


- 


- 


NGC 4309 


- 


- 




- 


NGC 4324 


150 


152 


169 


150 


NGC 4429 


276 


276 


294 


273 


NGC 4435 


237 


263 


266 


237 


NGC 4459 


244 


261 


273 


244 


NGC 4477 


249 


278 


295 


249 


NGC 4526 


328 


338 


348 


292 


NGC 4684 


136 


- 


137 


123 


NGC 4694 


96 




99 


93 


NGC 47 10 


175 


161 


202 


176 


NGC 4753 


301 


- 


302 


274 


NGC 5173 


180 




202 


166 


NGC 5273 


127 


- 


148 


125 


NGC 5379 


152 




158 


147 


NGC 6014 


148 


149 


185 


148 


PGC 29321 


154 




165 


154 


PGC 56772 


152 




153 


133 


PGC 581 14 










PGC 61468 


166 




166 


160 


UGC 05408 


139 




163 


133 


UGC 06176 


184 


196 


219 


184 


UGC 095 19 


160 


165 


182 


155 



Notes: This table contains the circular velocities derived from the JAM 
models, extracted at different radii. Galaxies that were removed from the 
ATLAS 3D sample have no JAM model available. Column 2 contains the 
circular velocities extracted at one effective radius. Column 3 contains the 
circular velocities extracted at the maximum radial extent of the CO, or 
12", whichever is smaller. Galaxies without interferometric observations 
are not included. Column 4 contains the maximum circular velocity. 
Column 5 contains the median of velocities after the peak in the circular 
velocity curve out to the maximum radius at which the model is constrained 
by the observational data. Some galaxies had insufficient data to constrain 
such a median between these two points, and are excluded. The error on 
these JAM velocities is estimated to be f»8%. This is a combination of an 
average 5% error in the observed stellar velocity dispersion constraining 
the models, and model errors of ss6% tCappellari et al. 2006). This does 
not include errors in inclination, which are harder to quantify. 




100 200 300 400 500 600 700 
w 2o,best/ sin (0 (km s" 1 ) 

Figure 4. Same as Figure[3] but for JAM model circular velocities measured 
at the maximum (top) and using the median velocity after the maximum 
(bottom). 



archival HI single-dish line-widths measured at 20% of the peak 
flux from HyperLEDA JPaturel et al.|2003| denoted W 20 ,hi). 

In Figure|5] we compare the single-dish HI line- widths to twice 
the projected JAM circular velocities measured at R c (top panel) 
and the observed CO line- widths (bottom panel). All galaxies from 
our sample with archival HI data at the time of retrieval are in- 
cluded. We find that although the average difference between the HI 
line-width and the JAM circular velocities is small (0.02 dex), the 
scatter is large (RMS scatter of ~0.1 1 dex), and the single-dish HI 
line- widths do not correlate statistically with the galaxy luminosity. 
The same is true when one compares with observed CO line- widths, 
with a larger scatter (~0. 15 dex). There is also a systematic trend in 
the residuals in the bottom panel of Fig. [5] where the galaxies with 
high CO rotation velocities, which are likely of high masses, have 
systematically smaller HI line-widths. 

Figure|5]suggests that some of the detected HI sources, at high 
CO velocities (and hence likely at high mass), are either not related 
to the galaxy under study or that the HI has been kinematically dis- 
turbed. This is consistent with the results of Morganti et al.|p006) , 
| William s et al. (20101 and Serra et al., (in preparation). 

Single-dish HI Line-widths in the LEDA database are how- 
ever known to suffer from many issues, including major problems 
with source confusion. Outer HI disks are also usually warped with 
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3.0 



2.0 2.2 2.4 2.6 2.8 3.0 

Log 10 (W 2asD km s" 1 ) 

Figure 5. Top: Difference between the deprojected HI line-width and twice 
the JAM model circular velocity measured at one effective radius, for galax- 
ies with archival HI data in LEDA, as a function of the model velocities. 
The solid line indicates where the JAM and HI velocities are identical. The 
dashed line is the mean difference, with the RMS scatter indicated in grey. 
Bottom: As for the top panel, but comparing observed CO single-dish and 
HI line-widths, for galaxies with boxy CO profiles only. 



respect to the inner regions, and hence a simple inclination esti- 
mates from optical images will be insufficient. Interferometric HI 
data allows one to identify relaxed disks, and in these rotation and 
the TFR can be studied (e.g. Morganti et al. 2006, Oosterloo et al. 
|2Tj07)|We1jrnans et al.|20"08"lloosterloo et al.|2010| >. To further in- 
vestigate the TFR, and the disk-halo consipracy in these galaxies, a 
future paper in this series will study the HI TFR of ETGs at large 
radii using superior interferometric HI data. 

4 Results 

Using the data presented in Table [T] we construct a series of 
Tully-Fisher relations to explore the effect of the different incli- 
nation and velocity measurements, denned in order of increasing 
complexity. The results are shown in Figure[6] The general form of 
the K B band Tully-Fisher relation we have adopted is 

M J , = a[log 10 ( i ^ T ) -2.6]+b, (2) 

where a is the slope and b is the zero-point of the relation. 

The package MPFIT < |Markw ardt 20091 was used to fit the in- 
verse of Equation [2] (regressing the observed rotation velocities) 
with the addition of an intrinsic scatter, which we iteratively ad- 
just to ensure a reduced \ 2 ~ 1- Full details of the fitting procedure 
can be found in |Williams eTal1 ( |20Trj) . 



Each plot in Figure [6] shows two best-fit lines. For the first 
we allow both the intercept and slope to vary, while for the sec- 
ond we fix the slope to that found by Tully & Pierce (2000]!, a = 
-8.78. For reference we also plot the A'-band spiral galaxy TFR of 
Tully & Pierce ( 2000) (see Equation[3j, which was constructed us- 
ing HI line-widths of spiral galaxies in 12 well separated clusters, 
and the i^-band SO TFR from |Williams et al.|p5~R)l l (see Equation 
|4j, which was derived from Jeans modeling of major-axis stellar 
kinematics from 14 edge-on SO galaxies: 

M K = -8.78 [log 10 ( j^fr) - 2.5] - 23.17, (3) 

Mks = -8.56 [log 10 - 2.4] - 24.02 . (4) 

4.1 CO 1\illy-Fisher relations 

4.1.1 Inclinations from galaxy axial ratios 

The CO TFR de-projected using i b / a is presented in the top two 
panels of Figure|6] The parameters of the best-fit relations are listed 
in Table[3] The top-left panel shows all galaxies, with open squares 
for those which do not show a boxy CO profile or where the pro- 
file shape is unclear, and solid triangles for those with clear boxy 
profiles. Only a loose correlation is observed, with clear outliers, 
almost all of which are galaxies without double horned profiles. 
A much tighter correlation is observed once one removes galax- 
ies with non-boxy profiles, as shown in the top-right panel. This 
is physically motivated, as we believe the gas in these galaxies 
does not reach beyond the peak of the rotation curve. Once the 
non-doubled-horned galaxies are removed both best-fit lines agree 
within the errors with the TFR of Williams et al. (2010), provid- 
ing strong a posteriori evidence that the CO in these galaxies does 
trace the circular velocity in a way that is consistent with dynamical 
models. The remaining outlier is NGC 7465 which we will discuss 
in more detail below. If this outlier is removed the best fit relation 
follows that of |Williams eTai~| < [20To] > even more closely. The in- 
trinsic scatter of the best fit relation with a free, and a constrained, 
slope is respectively 0.57 and 0.63 mag, with total RMS scatters of 
0.69 and 0.78 mag. 

4.1.2 Inclinations from dust axial ratios 

In a similar way we can construct a TFR using the inclinations de- 
rived from fitting ellipses to the dust distributions highlighted in 
unsharp-masked images. The CO TFR deprojected using id us t is 
presented in the middle-left panel of Figure [6] Once again, within 
the quoted errors, the best-fit relations are consistent with the result 
of |Williams et al.| ( |2010| l. The parameters of the best-fit relations are 
listed in Table [3] The intrinsic scatter of the best-fit relation with a 
free, and a constrained, slope is 0.65 and 0.78 mag respectively, but 
with a total scatter of 0.75 and 0.94 mag, respectively. The outlier 
is NGC 7465 - without this galaxies the observed intrinsic scatter in 
the unconstrained fit decreases to 0.46 mag, with a total scatter of 
0.6 mag, confirming the visual impression that the correlation using 
dust inclinations is tighter. Removing this galaxy also improves the 
agreement between our best-fit slope with that found by |Williams| 
|etall ( [20l0l >. 

4.1.3 Inclination from molecular gas modeling 

Using the inclinations from the tilted-ring models and the total ve- 
locity widths calculated from the interferometric data, we can con- 
struct yet another TFR. First we do this naively, using only the 
single-dish boxy classifications. This results in the middle-right 
panel of Figure [6] Some of the interferometric observations have 
insufficient spatial resolution to strongly constrain the inclination 
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of the molecular gas, and this is reflected in the error bars. As can 
clearly be seen in this panel, the best-fit relation is only marginally 
consistent with that of Williams et al. (20101. The intrinsic scatter of 
the best-fit relation with a free, and a constrained, slope is 0.43 and 
0.62 mag, with a total scatter of 0.60 and 0.82 mag, respectively. 

The interferometric data, at least in some cases, allow us to fur- 
ther refine our technique for identifying galaxies that have reached 
beyond the peak of the circular velocity curve, by identifying nu- 
clear rings and finding galaxies where a turnover can be seen in 
the position-velocity diagram. Interferometric observations are es- 
pecially powerful in this regard, as the extent of the CO can be 
compared to the expected circular velocity curve from dynamical 
models. For example, NGC7465, the obvious outlier in the pre- 
vious plots, features a misaligned circum-nuclear CO distribution, 
which has a double horned profile, but has no turnover in its veloc- 
ity field, and is extremely unlikely to reach beyond the peak of the 
galaxy rotation curve. 

Aperture synthesis observations also allow us to identify cases 
where a galaxy is misclassified using the single-dish data because a 
pointing error and/or extended molecular gas results in only some 
smaller portion of the total velocity width being detected. Such 
cases include NGC 4324, NGC 4477 and PGC 58 1 14. In a few cases 
where the profile shape was uncertain (due to low signal to noise 
single dish spectra) the interferometric data reveals that the gas does 
extend beyond the peak of the rotation curve, allowing us to include 
them. This is the case for NGC 3489, NGC 4292 and UGC 095 19. 

The TFR in the middle-right panel of Figure [6] can thus be 
improved using the full knowledge gained from the interferomet- 
ric data (discussed above), and this results in the TFR shown in the 
bottom-left panel of Figure|6] The best-fit TFR using this additional 
information agrees with that of Williams et al. (2010). The parame- 
ters of the best-fit relations are again listed in Table [3] The intrinsic 
scatter of the best-fit relation with a free, and a constrained, slope 
is 0.32 and 0.33 mag, with total a scatter of 0.54 and 0.58 mag, 
respectively. 



4.1.4 Combining interferometric and single-dish data 

Finally, we can combine our single-dish and interferometric obser- 
vations, using the hybrid inclinations and velocities described in 
Section l2.6.1l These choices lead to the TFR shown in the bottom- 
right panel of Figure [6] The parameters of the best-fit TF relations 
are liste d in Table [3] Th e best-fit relations closely agree with the 
result of Williams et al. (2010), and the fits with a free and a con- 
strained slope have an intrinsic scatter of 0.36 and 0.37 mag, with a 
total scatter of 0.54 and 0.57 mag, respectively. The intrinsic scat- 
ter found in this relation is similar to that found by a large study 
of the intrinsic scatter of the spiral galaxy TFR (<7i nt ~ 0.4 mag at 
bands between g and z; |Pizagno et al. 2007[>, and only very slightly 
larger than that found by Williams et al. 1 [2010[ > (er^t ~ 0.3 mag 
at A's-band). The small total scatter reflects both the better number 
statistics and the use of the best possible measures of rotation and 
inclination. 

Using this hybrid relation we can investigate if the residu- 
als have any systematic behavior with respect to the properties of 
the hosts. We find that the residuals are statistically uncorrelated 
with the gal axy if B -band luminosity jSkrutskie et al.|2006) , optical 
colour (ff-r; |Adelman-McCarthy et al.|2008 Scott et al., in prepara- 
tion), velocity dispersion, galaxy total dynamical mass (Cappellari 
et al. in preparation) and spin parameter \r (Emselle m et al.|2010| 
Paper III) 



5 Discussion 

5.1 CO as a tracer of galactic potentials 

The results presented in Section [3~T] show that CO emission accu- 
rately traces the potential of early-type galaxies, if one has a suit- 
able method of identifying which galaxies have extended molecu- 
lar discs. Our work demonstrates that selecting galaxies via profile 
shape is robust and does not systematically affect the measured ro- 
tation velocity. The main possible contaminant introduced by this 
selection is galaxies that have molecular rings, which would display 
a boxy profile without necessarily rotating at Vn at . Circumnuclear 
rings pose the greatest problem, as they are very unlikely to rotate 
at Vfl a t, unlike inner and outer rings which are located at resonance 
points within the flat part of the galaxy circular velocity curve (e.g. 
|Buta|19 86). The number of galaxies with circumnuclear rings in our 
sample is likely to be small, as their presence in significant numbers 
would bias the best-fit line in Figure[5]to small values. This conclu- 
sion is also supported by the interferometric survey of Alatalo et 
al., (in preparation), which suggests that only a small percentage of 
local early-type galaxies have their molecular gas restricted to the 
circumnuclear region. 

One should also remember that the circular velocity curves of 
ETGs rise very quickly, normally peaking within one effective ra- 
dius (e.g. Williams et al. 2009). This allows the molecular distribu- 
tion to be compact, but still have enough gas beyond the peak of the 
galaxy circular velocity curve so that the linewidth is a good proxy 
of the circular velocity in that region. The good agreement between 
our derived TFR and that of Williams et al. (2010) suggests that the 
disk-halo conspiracy holds for our early-type galaxies, at least out 
to 3R C . If this conspiracy holds even further out into the dark mat- 
ter dominated parts of these galaxies (which some authors suggest 
it does not, e.g. Dutton et al. 2010) will be considered in a future 
paper in this series. 

5.2 Inclination measures 

The problem of measuring the inclination of early-type galaxies is 
not new, but it is complicated here as the CO is often not aligned 
with the sta rs (e.g.|You ng 2002 ; Schinnerer & Scoville 2002 ; Young| 
|et al.|2008| [Crocker et al.|2008||2009||2010> , and it is the CO in- 
clination that we require. As we have shown, methods using galaxy 
optical axis ratios are likely to yield large uncertainties and increase 
the scatter in the TFR, but they do not seem to significantly affect 
the determination of its slope or zero-point. This is likely because 
errors in inclination are most problematic when the galaxy is close 
to face-on, and the number of nearly face-on galaxies in our sample 
is small. The profile shape selection criterion also actively selects 
against face-on galaxies, where the velocity width is less than a few 
channel widths (no matter how extended the molecular gas is). In 
many HI TFR papers, galaxies with an inclination less than 45° 
are discarded (e.g. Tully & Fisher|1977) , so this is not especially 
worrying. If we apply this cut in inclination to our sample, we re- 
tain within the errors the same Tully-Fisher slope and zero-point, 
and the RMS scatter decreases. Even with a very harsh cut of 65° 
we retain the same relation. This demonstrates that our results are 
robust, and the slope and zero-point are minimally affected by the 
choice of inclination correction. 

As hinted above, while the slope and zero-point are robust, 
the scatter in the relation is affected by the inclination estimation 
method chosen. Moving from galaxy axial ratios to inclinations de- 
rived from dust features decreases the RMS scatter around the best- 
fit relation by ~0.15 mag, if one removes the one obvious outlier. 
Our results do reveal a reduction of the total scatter when the in- 
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Figure 6. CO TFRs constructed using various inclination and line-width measurements. The correlations shown in each plot and indicated in the legend are the 
spiral galaxy TFR of Tully & Pierce (2000 1, derived from HI line-widths (short dashes), and the SO TFR of Williams et al. (2010 1, derived from Jeans modeling 
(long dashes). Also plotted are the best-fit relations for our data, of the form given in Equation[2] for both an unconstrained fit (dot-dashed) and a fit where 
the slope is fixed to that found by Tully & Pierce (20001 (triple dot-dashed). The best-fit parameters are listed in Table[5] Top-left: Inclinations derived from 
galaxy optical axial ratios (see Section [2.5.1) . CO detections with boxy profiles are shown as filled diamonds, galaxies with non-boxy profiles and uncertain 
classifications are shown as open squares. Top-right: As top-left, but galaxies with non-boxy profiles are removed. The outlier is NGC7465. Middle-left: As 
top-right, but using inclinations from unsharp-masked dust images (see Section |2.5.2[ . The outlier is again NGC7465. Middle-right: As top-right, but using 
interferometric W20 values and inclinations derived from tilted-ring fits to the interferometric data (see Section |2.5.3) . Bottom-left: As middle-right, but using 
the interferometric data to identify galaxies with CO which reaches beyond the peak of the galaxy rotation curve. Bottom-right: A hybrid plot using the best 
inclinations and velocity widths as defined in Section p.6. 1| Both best-fit lines overlay the relation of Williams et al. 



clinations are determined from fitting a thin disc model to interfer- 
ometric data, despite our limited number statistics. With a larger 
sample one might expect the scatter to decrease further. In both 
cases the decrease in scatter is likely due to properly accounting 
for molecular gas distributions that do not rotate in the same plane 
as the stars. Unsurprisingly, to obtain the tightest relation one must 
use measures of inclination (such as unsharp-masked dust images or 
interferometric imaging of the molecular gas itself) that truly trace 
the inclination of the gaseous component. 

5.3 Intrinsic scatter of the CO TFR 

The intrinsic scatter of our hybrid TFR is 0.36 mag, broadly con- 
sistent with typical values found in dedicated studies of the spiral 
galaxy TFR in the optical, such as those by Pizagno et al. (2007 ) 
(0.4 mag at g- and 2 -band) and |Kannappan et~ (2002) (0.4 mag 
at r-band). Our intrinsic scatter is slightly larger than that found 



for SO galaxies by |Williams et al.| ( |2010l > («0.3 mag at K s -baa&). 
One reason for the increased scatter may be that Willia ms et al.| 
( [20 1 0| l considered SO galaxies only, whereas our sample includes 
both morphologically classified elliptical and SO galaxies, that span 
a wider range in Ar, the stellar specific angular momentum (Em- 
|sellem et al.|2007) . This increased scatter may therefore be due to 
spanning an increased range in internal dynamics. 

We do not have sufficient number statistics to constrain an el- 
liptical galaxy TFR separately, but the earliest type galaxies (se- 
lected by optical morphology) do lie on average to the faster and/or 
dimmer side of our relations. There is however no systematic be- 
havior within the residuals to the fit as a function of Ar. 

5.4 Offset from the spiral TFR 

Exploring changes in the M/L of galaxies as a function of morpho- 
logical type is one of the main motivations for studying the early- 
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Table 3. Best-fit parameters of the TFRs shown in Figure[6] 


Inc. method 


a (mag) 


b (mag) 


cr int (mag) 


"total (mag) 


W 2 0,sd 


W20,inter 


Selection 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Stellar bla 


-5.61 ± 0.71 


-23.62 ± 0.16 


0.75 


0.82 


X 




none 


Stellar bla 


- 


-24.01 ± 0.23 


1.30 


1.42 


X 




none 


Stellar bla 


-7.82 ± 1.12 


-23.27 ± 0.15 


0.60 


0.69 


X 




SD 


Stellar bla 




-23.25 ± 0.16 


0.66 


0.78 


X 




SD 


Dust 


H If, _1_ 1 1 O 

- /.jo ± i.iy 


-LJ.Dl ± U.iO 


U.oo 


u. / o 


X 




en 

aij 


Dust 




-23.30 ±0.19 


0.79 


0.94 


X 




SD 


Tilted-ring 


-6.68 ± 1.21 


-23.26 ±0.17 


0.48 


0.63 




X 


SD 


Tilted-ring 




-23.23 ± 0.23 


0.65 


0.86 




X 


SD 


Tilted-ring 


-8.07 ± 1.10 


-23.05 ± 0.12 


0.32 


0.54 




X 


Inter 


Tilted-ring 




-23.03 ±0.13 


0.33 


0.58 




X 


Inter 


Hybrid 


-8.38 ± 0.70 


-23.12 ±0.09 


0.36 


0.54 


X 


X 


Inter+SD 


Hybrid 




-23.11 ± 0.09 


0.37 


0.57 


X 


X 


Inter+SD 



Notes: Column 1 describes the inclination method used to deproject the measured line-width. Columns 2 and 3 show the best-fit slope and 
zero-point, respectively, of the Tully-Fisher relation of the form Mjc s = a[log 10 (W20 /km s -1 ) — 2.6] ±b. Where Column 2 is blank the fit 



was constrained to have the same slope as that found by Tully & Pierce 1 2000 1, a=-8.78. Column 4 shows the intrinsic scatter (<r int ) required 
by the fit to produce a reduced x 2 ~ 1- Column 5 shows the total RMS scatter around the best-fit relation. A cross in Column 6 denotes that the 
velocity widths used are from single-dish data, while a cross in Column 7 denotes that interferometric line- widths were used. Column 8 lists the 
selection method used to find galaxies that have gas beyond the peak of the circular velocity curve. 'SD' means the galaxies single-dish profiles 
are boxy, while 'Inter' means that interferometric data were used to determine if the gas reaches beyond the peak of the circular velocity curve. 
The hybrid relations use both single-dish and interferometric data, as described in Section [2.6.1| 



type TFR. Willi ams et al.| (|2010) have highlighted the importance 
of comparing TFRs for each morphological type derived using the 
same tracer and methods, in order to minimize systematic differ- 
ences. In this spirit we use the spiral galaxy line-widths published 
by |Chung et al.|(2009| >, measured from on-the-fly mapping of 18 
Virgo Cluster members, in order to define a spiral CO TFR. 

We construct a Tully-Fisher relation for the Chung ^t"aL] ( |2009} 
spirals in the same way as for our early-type sample, using 2MASS 
K B> total magnitudes and their W 2 o line-widths corrected for incli- 
nation. These inclinations were derived from the galaxy axial ratio 
and the morphological type using the classical Hubble formula from 
|Hubble| ( |l926) . This formula effectively differs from Equation |2.5.1| 
by making qo a function of morphological type. For spiral galaxies 
such as those in the |Chung et al.| ( f2009^ sample this correction is 
small, corresponding to a qo less than 0.05, and hence this differ- 
ence is unlikely to systematically effect the comparison of our re- 
sults. We differ from Chung et al. ( 2009 ) however by not assuming 
a set distance to the Virgo Cluster, which has significant substruc- 
ture. Instead we use the distance to each galaxy as determined for 
the ATLAS 3D parent sample (Paper I). The resulting TFR is shown 
in Figure [7] The values used to create this figure are listed in Table 

s 

The best-fit spiral TFR with an unconstrained gradient is 

W 20 



M Ks = (-10.8±2.5) [lo gl 



/ Wgo \ 
Vkms- 1 / 



-2.6]-(24.1±0.2) ,(5) 



where the intrinsic scatter about the correlation is 0.68 mag, with a 
total scatter of 0.72 mag. The TFR obtained by constraining the fit 
to have the same slope as that of Tully & Pierce ( 2000 I is 



Mi 



8.78 [log 10 (^pr) - 2.6] - (24.02 ± 0.14) . 



(6) 



with an intrinsic scatter of 0.55 mag and a total scatter of 0.59 mag. 

Within observational errors, the best-fit free and constrained 
spiral CO TFRs are consistent with the relation of Tully & Pierce 
(2000). The offset between the unconstrained Chung et al. (2009) 
spiral CO TFR and the hybrid early-type CO TFR (discussed in 
Section |4.1.4[ ) is 0.98 ± 0.22 mag, consistent within errors with 
the value found by Bedregal et al. (2006) and the offset found by 



Table 4. List of the parameters used to create the spiral TFR in Figure|7] 



Galaxy Type 



(1) 



(2) 



Distance 

(Mpc) 

(3) 



niKs 
(mag) 
(4) 



(deg) 
(5) 



W20 
(kms" 1 ) 
(6) 



NGC4254 


Sc 


16 


5 


6 


93 


29 


± 


7.3 


221 


± 


1 


NGC4298 


Sc 


16 


5 


8 


47 


59 


± 


2.4 


270 


± 


4 


NGC4302 


Sc 


14 




7 


83 


90 


± 


0.1 


353 


± 


1 


NGC4303 


Sc 


16 


5 


6 


84 


19 


± 


12. 


162 


± 


2 


NGC4321 


Sc 


15 


85 


6 


59 


38 


± 


5.2 


239 


± 


2 


NGC4402 


Sc 


16 


5 


8 


49 


80 


± 


0.7 


267 


± 


3 


NGC4419 


Sa 


13 


12 


7 


74 


82 


± 


0.6 


318 


± 


2 


NGC4438 


Sb 


16 


5 


7 


27 


87 


± 


0.2 


272 


± 


5 


NGC4501 


Sbc 


15 


3 


6 


27 


60 


± 


2.3 


518 


± 


2 


NGC4527 


Sb 


14 


05 


6 


93 


75 


± 


1.1 


376 


± 


2 


NGC4535 


Sc 


15 


75 


7 


38 


41 


± 


4.7 


252 


± 


5 


NGC4536 


Sc 


14 


09 


7 


52 


59 


± 


2.4 


318 


± 


2 


NGC4548 


Sb 


18 


71 


7 


12 


35 


± 


5.8 


234 


± 


6 


NGC4569 


Sab 


16 


5 


6 


58 


69 


± 


1.6 


311 


± 


1 


NGC4579 


Sab 


21 


3 


6 


49 


39 


± 


5.0 


329 


± 


1 


NGC4647 


Sc 


16 


5 


8 


05 


34 


± 


6.0 


156 


± 


1 


NGC4654 


Sc 


16 


5 


7 


74 


58 


± 


2.5 


285 


± 


1 


NGC4689 


Sc 


16 


5 


7 


96 


39 


± 


5.0 


182 


± 


1 



Notes: The galaxy sample listed in Column 1 comes from Chu ng et al.| 
{2009). Column 2 lists the spira l galaxy type following the Hubble scheme, 
(RC3;|de Vaucouleurs et al. 19911. The distance to the galaxy in Column 3 
is from the ATLAS 30 parent sample (Paper I). These distances are drawn 
preferentially from |Meiet "51^2007) and |Tonry etal~|j2001) . Column 4 lists 
the apparent Kg -band magnitudes, taken from 2MASS. Column 5 is the 
optical galaxy axis ratio, taken from NED, as described in Section [2.5. 1| 
The line-widths in Column 6 are taken directly from Chung et al. ( 2009). 



Williams et aT]p0"T0) between their SO TFR and the |Tully & Pierce| 
(2000) relation. 

This result is nevertheless surprising since Will iams et al.| 



( |2010| l argue that their offset with |Tully & Pierce|p000| > is flawed 
by systematic effects (the tracers and methods being different), and 
they prefer the much smaller offset (0.5 ± 0.14 mag) from their 
spiral galaxy sample treated identically to the SOs. 

It is not immediately clear how our own result can be recon- 
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Figure 7. Spiral galaxy CO TFR built in the same way as that shown in 
Fig. [6] but using W20 values and inclinations from the work of Chu ng et al.| 
p09) . 



ciled with this, as we use identical tracers and methods as |Chung| 
|et al. | {20091. However, we suspect that the different results are due 
to the mix of morphological types in the spiral samples. |Williams] 
et al. (2010) use a spiral sample consisting mostly of Sa-Sb spi- 
rals, while the sample of Chung et al. ( 2009 1 includes a larger frac- 
tion of later-type galaxies (mostly Sc), which are known to have 
a larger zer o-point (e.g. | Roberts| 1978} |Rubin et al.|1985}|Masters| 
|et al.|2008[|Shen et al.|2009| >. Separating the |Chung et al.| < |2009) 
sample more finely into morphological types does not necessarily 
support this explanation, but the number statistics are very limited. 

5.5 Origin of the offset TFR for ETGs 

In the ETG sample used in this paper, a few galaxies are consistent 
with the spiral galaxy relation and many lie between the spiral and 
early-type relations, consistent with the predictions of passive fad- 
ing models (e.g. |Dressler|1980| [Dressier et al.|1997| >. The fact that 
molecular gas is present in all these galaxies, however, suggests that 
many must have ongoing star formation (even if weak), and hence 
evolutionary models relying on a purely passive evolution are ruled 
out. 

The presence of molecular gas mass fractions ranging from 
<0.1 to 5% in these galaxies perhaps suggests an answer to some of 
the problems with passive evolution models. Bedregal et al. (2006) 
and Willia ms et al.| (2010) both discuss the timescale implied by 
the ETG-spiral TFR offsets observed, assuming passive fading, and 
find it uncomfortably short and inconsistent between the B and K s 
bands. The presence of low levels of residual star formation would 
prolong the period it takes for galaxies to fade, and would contribute 
relatively more light at blue wavelengths than in the infrared. The 
molecular gas mass fraction in these galaxies is however small , so 
it is unclear if such a low mass fraction of young stars could extend 
the fading timescale sufficiently. 

An alternative explanation of the TFR offset is that the mor- 
phological (and luminosity) evolution is accompanied by a change 
in the size of the galaxies, affecting the galaxy rotation rather than 
its luminosity. Williams et al. (2010) suggest that a systematic con- 
traction as spirals transform into SOs, consistent with the trend with 
morphological type of the size-luminosity relation ( Courte au et al.| 
2007 1, could explain the offset from the spiral TFR. Clearly, a com- 
bination of both effects is also possible. 



5.6 High-mass end of the CO TFR 

As discussed in the introduction, one motivation for studying the 
CO TFR is the ability to probe the full mass range of galaxies, 
from small disc galaxies to large spheriods, with the same tracer 
and hence with the same assumptions and systematics. As an ex- 
ample of the power of this approach, we can investigate if our sam- 
ple shows any evidence of a change in the slope of the TFR at the 
high-mass end, as repo rted by some authors (e.g. [P eletier & Willner| 
|T993]|Verheijen|2001||Noordermeer & Verheijen|2007) . 

Our sample contains eight galaxies above the luminosity 
threshold suggested in previous works (Mx s ^-23.75). As can be 
seen in the bottom-right panel of Figure[6] all but one of these sys- 
tems are on, or to the right of, the best-fit relation for the whole 
sample. When one compares the best-fit TFR derived for the whole 
sample to that for just the fainter galaxies (Mk s ^-23.75), there 
is a small difference in the gradient and zero-point, in the direc- 
tion found by other authors. We do not have the number statistics 
to properly investigate the size of this effect, however, and hence 
only conclude that our results are consistent with the presence of a 
different slope at the high-mass end of the TFR. 

5.7 Single-dish vs aperture synthesis 

We have demonstrated in this work that single-dish spectra are suffi- 
cient to construct a robust CO Tully-Fisher relation, and that galaxy 
inclinations can be taken from optical imaging without noticeably 
biasing its slope or zero-point (Section [7!2| . However, one must ex- 
ercise caution, as the small beams of millimeter telescopes can oc- 
casionally underestimate the line-widths of galaxies with extended 
molecular discs. Barring beam dilution effects, one would expect 
single-dish profiles to perform even better at higher redshifts, where 
the angular size of the target galaxies decreases. One caveat to this 
is that as the angular size of galaxies decreases, source confusion 
may become a problem. 

Naively, interferometric observations should be the best way 
to obtain an estimate of the true inclination of the molecular gas, 
and hence reduce the artificial scatter in the TFR. Aperture syn- 
thesis observations also tend to have a large field of view, and are 
thus less likely to face problems with extended gas distributions. 
Interferometric observations do, however, usually require a larger 
investment of observing time per source. One must also choose the 
angular resolution of the data carefully, in order to avoid resolving 
out the most extended structures, while still having sufficient spa- 
tial information to constrain a tilted-ring model. Both single-dish 
and interferometric studies of the TFR should thus have a bright 
future. 

6 Conclusions and Future Prospects 

In this paper we have presented the first CO Tully-Fisher re- 
lation for early-type galaxies. We have shown that CO line widths 
are a reliable tracer of the circular velocity in fast-rotating ETGs, 
and relations derived in this way agree well with ETG TFRs de- 
rived by other authors. We showed that our CO TFR relation has 
a robust slope and zero point, and a comparable scatter to that of 
TFRs derived with other techniques. The power of the CO Tully- 
Fisher relation is the ability to perform the same analysis easily for 
both spiral and early-type galaxies (i.e. identical tracer and method), 
with purely observable quantities. This technique is therefore par- 
ticularly suitable for comparing the morphological variations of the 
zero-point of the Tully-Fisher relation. 

The results presented in this paper show that early-type galax- 
ies follow a TFR that is offset by nearly one magnitude at i^ s -band 
from that of spiral galaxies. This offset is similar to that found 
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by previous authors. An offset ETG TFR is usually thought to be 
caused by passive evolution of spiral galaxies into lenticulars, af- 
ter an abrupt cessation of star formation. However, the presence of 
molecular gas in these ETGs suggests that the majority of this sam- 
ple must have ongoing residual star-formation, such as that revealed 
in a subsample of these galaxies (e.g. |Temi et a l. 2009 ; Shapiro 
|et al.|2010||Crocker et al.|2010t . This should increase the timescale 
to fade by the required ~1 mag. Some models, such as those that 
involve a change in the size of galaxies as they transform, do not 
require star formation to completely cease, providing an alternative 
explanation that is fully consistent with our data. Unfortunately the 
low molecular gas fractions in our systems make drawing any firm 
conclusions difficult. 

Another area where the CO TFR (for all spectral types) ex- 
cels is the ability to extend the technique to higher redshifts, as 
discussed in Tutui et al. (2001 ). HI observations are currently lim- 
ited to z < 0.2 (e.g.|catinella|2008|l, and stellar kinematics to z < 2 
(e.g. |van Dokkum et al.||2009| |Cappellari et al.||2009[ |Cenarro &| 
Trujillo 2009} whereas CO has already been detected in galaxies at 
redshifts z > 6 (e.g. |Walter et al.|2004| >. The molecular content of 
galaxies is also expected to increase with redshift, while the HI con- 
tent stays approximately constant (Obrcschkow & Rawlings 2009, 
|Bauermeister et al.|2009"l |Daddi et al.|2010[|Tacconi et al.|2010) . 
CO velocity fields can be resolved, and the rotation (and hence the 
dynamical mass) estimated even out to high redshift (e.g. Walter 
et al. 2004; Carilli et al. 2010). This raises the exciting prospect 
of being able to directly track the mass-to-light ratio evolution of 
galaxies as a function of both redshift and mophological type, us- 
ing an easily observable tracer (and without requiring dynamical 
modeling). New optical/IR instrumentation is just starting to allow 
surveys of galaxy stellar kinematics out to z ~ 2 (e.g. Cappellari 
et al. 2009; Forster-Schreib er et al.|2009) , and comparison of these 
samples to line-widths derived with CO is a promising immediate 
avenue for future studies of the TFR. 

|Tutui et aT7| < |200l| l showed that constructing a (spiral galaxy) 
CO TFR is already possible at up to redshifts z~0.1, but with the 
next generation of mm-wave telescopes due to enter service over 
the next few years, we should gain access to CO lines across the 
redshifted universe for all morphological types. The Redshift Re- 
ceiver System (RRS) on the Large Millimeter Telescope (LMT), 
for instance, is optimized for the detection of redshifted transitions 
of the CO ladder from star-forming galaxies at cosmological dis- 
tances, and will generate a large database of CO line-widths. If one 
is able to optically image these galaxies to estimate their inclination 
and morphology, then this will become an excellent resource for 
Tully-Fisher analyses. With the ability to detect spectral line emis- 
sion from CO in a galaxy like the Milky Way at a z~3 in less than 
24 hours, ALMA will also be in an excellent position to study the 
CO TFR at high redshift. The higher frequency bands that will be 
available at ALMA also open up the possibility of detecting [CII] 
emission (rest wavelength 158 /im), which is brighter than CO in 
many galaxies in the local universe < |Crawford et al.|1985|> an d pos- 
sibly enhanced further at high redshifts jMaiolino et al. ;2009 1. [CII] 
may therefore become an important tracer of the circular velocity of 
high-redshift systems. 

Although technical aspects will certainly be improved and the 
interpretation of this relation remains uncertain, we have shown 
that the ETG CO TFR is a tantalising and promising tool for galaxy 
evolution. 
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